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Abstract 

We discuss the evaluation of Luria-Delbriick fluctuation experiments un- 
der Bellman-Harris models of cell proliferation. It is shown that under 
certain very natural assumptions concerning the life-time distributions 
and the offspring distributions of mutant and non-mutant bacteria, the 
suitably normed and centered number of mutants contained in a large 
culture of bacteria converges to a stable random variable with index 1. 
The result obtains under the assumption that the mutation under consid- 
eration is 'neutral' in the sense that on average, mutant cells produce the 
same number of offspring as non-mutant cells. Thus, provided this condi- 
tion is met, fluctuation experiments can be evaluated without knowledge 
of how cell proliferation proceeds in detail. This resolves a long-standing 
issue in the theory of Luria-Delbriick fluctuation analysis. 
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newal Theory 

Mathematics Subject Classification (2000): 60J85, 92D15 



1 Introduction 

Fluctuation analysis, or the art of measuring mutation rates by means of Luria 
and Delbriick's fluctuation test, has come into prominence sixty-odd years ago 
P|, when it was realized that the huge fluctuations of the number of phage- 
resistant bacteria in sensitive cultures of Escherichia coli could be taken as 
evidence that the phage-resistant bacteria had arisen in these cultures for no 
particular reason: For if they had, what better reason could there be for a 
bacterium to become phage-resistant than the encounter with the phage itself? 
But then comparable cultures of bacteria should, upon exposure to comparable 
amounts of the phage, produce more or less the same number of mutant colonies 
on solid medium. They don't, even if one allows for some variability in the 
number of mutant colonies a given culture of bacteria produces in an actual 
experiment, and so the development of phage- resistance by these bacteria is not 
triggered by the phage. Indeed, under the assumption that resistant bacteria 
are only produced in the presence of the phage, the actual number of mutant 
colonies emerging from reasonably many cultures of bacteria on solid medium 
should be distributed over the cultures in a Poisson-like manner. This is not the 
case, and so there must be at work some other mechanism by which sensitive 
bacteria are turned into resistant ones. Fluctuation analysis aims to illuminate 
this mechanism by way of a thorough analysis of the distribution of mutant 
colony count in what has come to be called a fluctuation experiment, which is a 
fancy term for the fact that a number of bacterial cultures have been grown from 
an initially small number of non-mutant bacteria, and the number of mutants 
each culture contains has been recorded (see, e.g., [21001E1II300III1 for 
more or less recent reviews of how fluctuation experiments are performed in 
the laboratory) . A particular emphasis of this analysis lies on determining the 
probability p (not a standard notation) for a cell to undergo mutation during 
division. It is customary to refer to this probability as the mutation rate. 

It is true that the fluctuation-analytic approach has faced recent competi- 
tion from certain other, experimentally more sophisticated methods jlOllllllT^l 
II HI 1141 IT5] . Still, it may serve as a welcome complement to these, all the more 
as certain parameters of possible interest besides the mutation rate can be de- 
termined from a fluctuation experiment as well |16|. For this, one requires a 
mathematical model of how to calculate the so-called Luria- Delbriick distribu- 
tion of mutant cells an individual culture of bacteria contains just before it is 
plated on solid medium. Now, it is not only since the (relatively recent) induc- 
tion of fluctuation analysis into the 'molecular biology wing of the Museum of 
Elegant Science' (the copyright for this nice one goes to Patricia Foster |12j) 
that questions concerning the mathematical theory of bacterial populations be- 
coming enriched with mutants have been paid a certain amount of attention. 
Not a very big lot of it, mind you, but the occasional article on the theory of 
the Luria-Delbriick distribution has popped up in the journals for longer than 
the recent half of a century (e.g., [T71[iTn[TT?10r)linil^ll^ll24ll^ll^)ll77l[^] 
besides those already quoted), and will, in all likelihood, continue to do so un- 
til the community is agreed that the issue simply isn't interesting anymore. 



To the author's knowledge, however, the influence that different models of cell 
proliferation might exert on the calculation of the Luria-Delbriick distribution 
and the evaluation of a fluctuation experiment has never been touched upon 
in detail. What one needs to look at in particular is non- exponential life-time 
distributions of cells. This is what makes clonal expansion (if only cell numbers 
are considered and not, say, the distribution of cells over possibly a contin- 
uum of age-classes) non-Markov and, therefore, a little messy. True enough, 
the probability generating function (PGF) of the Luria-Delbriick distribution 
is, the non-Markovian character of bacterial proliferation notwithstanding, easy 
to write down. But what's next? Even with the help of a computer program 
for symbolic calculation, the actual computation of the PGF in question isn't 
all that easy. (The reader who disagrees is cordially invited to try her hand 
at Equation Q below - and no, it is not permissible to assume that life-time 
distributions are exponential!) Is it all worth it? 

It could be argued that, to be on the safe side, one could always conduct 
a fluctuation experiment such that a reasonable fraction po of cultures will 
not contain any mutants at all. Since the probability that an initially very 
small culture of bacteria will not have produced any mutants by the time it has 
reached a certain size n (which is to say, after ~ n cell divisions) is (1 — p) n , 
an obvious estimate for the mutation rate would be p ~ — logp^/n. This would 
be true regardless of how the proliferation of cells proceeds in detail. The 
performance of this so-called po-method, however, is generally poor [T^l, so that 
the question is whether other methods of evaluating a fluctuation experiment 
need to be adapted to non-Markovian models of cell proliferation. We will see 
that even for 'neutral' models of mutation, the answer is in general yes. Just 
how much these methods need to be adapted is the main result of this paper. A 
translation of what should be understood by 'neutrality' in the context of single 
cell populations is as follows: a mutation is 'neutral' if, on average and in the 
long run, a mutant cell produces as many offspring as a non-mutant one. We 
shall be more precise about this in the following section. 

2 The Luria-Delbriick Distribution 

We shall take it for granted that the PGF of the Luria-Delbriick distribution is 
of the form 



9ld(s) :=: ^2p r s r = exp (m(g(s) - l)j , (1) 

r=0 

where m is parameter that we shall soon return to, 



and 



g(s):=P e-^F u (s)du, (2) 

Jo 

F u (s) = (1 - G'(u))s + f f o F u _ y (s) dG'(y) (3) 

Jo 



is the PGF of particle numbers in a Bellman-Harris process with first-generation 
offspring PGF 



fc=0 



/( S ):=y> feS fc , (4) 



and life-time distribution function G'(-) I2H1- I n words, G'(i) is the probability 
that a newborn mutant cell lives at least until time t, and -Kk is the probability 
that once it splits into progeny, it will split into exactly k of these. We will 
always assume that proliferation is supercritical, that is, 1 < /'(l) =: p' < oo. 
This implies, in particular, that there exists q 6 [0, 1) such that f(q) = q. 
For practical purposes, of course, /(s) can be taken to be only a quadratic 
polynomial in s, but this doesn't make anything about the calculations easier. 
We shall give little justification for Equation (JIJ in this paper; it is in agreement 
with (and in fact a generalization of) previous work on the subject (notably 
|19l 1251 EH! 1321 \'6'6\ ) . and in any case seems reasonable enough: the point is that 
with the mutation rate p being generally small, the first mutants will arise only 
when the population as a whole is already well into the stage of Malthusian 
growth, so that the probability for a mutant to arise u units of time before 
the present moment is ~ e _/3 " smaller. Here and in Equation J5J), (3 is the 
Malthusian parameter for the growth of a population of exclusively non-mutant 
cells, which for our purposes is a number that satisfies 



/<■ 



e-^dG°(t) = l, (5) 



o 



where p° is the expected number of non-mutant first-generation offspring of a 
non-mutant cell, and G°(-) is the life-time distribution of non-mutant bacteria. 
We now define a mutation to be 'neutral' if f3 also satisfies 

/"OO 

H a / c-"*dG*(«) = l, (6) 

Jo 

where p* is the expected number of mutant first-generation offspring of a newly 
mutated bacterium. We assume that both \x* (we frequently write p instead 
of p') and fi° are finite, that neither G°(-) nor G'(-) are lattice-like, and that 
G°(0) = G°(Q) = 0. Then at each given moment in time, there are only finitely 
many bacteria (or 'particles') alive jHOI, and provided that the population has 
not become extinct (which is the case with probability 1 — q > 0, and which 
we may easily assume), the population grows essentially as ~ e^* with t. The 
expected number of mutations in the population grows as ~ pe' 3 *. One of the 
basic assumptions of fluctuation analysis is that, as one allows the population to 
grow for longer and longer periods of time, the mutation rate in the population 
has always been so small that 

pe 131 — > m 

as t — > oo, where m < oo. This is the parameter entering into Equation (0), 
and the one of predominate interest to fluctuation analysts. 



The problem is thus to extract m from experimental data under the model 
assumption J2J). To this end, we impose the following restriction on f(s): Define 

Ks) := IfM , (7 ) 

1 — s 

where /(s) is the PGF in J3J. Our standing assumption is that, for s in the 
neighbourhood of zero, 

/*• - h{\ - s) = s» (- log *)-°£(- logs) , (*) 

where either w > and a 6 (—00,00), or w = and a > 2. £(•) is slowly 
varying at infinity, and in case u> = and a = 2 must be such that the integral 

&* M 

is finite. (*) is slightly stronger than the 'X log X '-condition of the Kesten- 

Stigum theorem. Via Tauberian arguments, it translates into X)fc=« ^Vk ~ 

n u '~ 1 (logn)~ a C(logn) for w < 1, else there is little intuitive about it. It is 

inspired by a paper of Uchiyama 34 . 

We shall show that if (*) holds and if the mutation under consideration is 

neutral in the sense that the Malthusian parameter (5 is determined by either 

Equation (JSJ) or ©, the law in QJ is in the normal domain of attraction of a 

stable distribution with index 1. Now the model considered by Lea and Coulson 

|19| in their celebrated formulation of the LDD is one of a neutral mutation and 

meets (■*•) with h(s) = 1 + s. Therefore its law is attracted to a stable distribution 

with index 1, and any part of the statistical machinery developed in Lea and 

Coulson's paper which rests on this fact alone can be can be applied no matter 

how cell proliferation proceeds in detail. The catch is that one will, in general, 

not obtain an estimate for the mutation rate itself, but only for the mutation 

rate times a factor 

• /i — 1 , . 

Ul := Wffte-VdG'V) ' ( ' 

It is seen that n' is completely determined by the life-time distribution G'(t) 
and the expected number of offspring (/i = /!*) of mutant bacteria. If either is 
unknown, n* must be determined from auxiliary experiments. This may seem 
a drawback of the method advocated in this paper, but it is actually a problem 
inherent to fluctuation analysis itself: It seems intuitive that by mere inspection 
of the distribution of mutant colonies in an experiment, it might be difficult to 
disentangle how mutation and cell proliferation work together to produce the 
actual total number of mutants each population contains (think of a population 
where mutants are produced at a higher rate but it takes longer for mutant 
cells to divide). From this point of view, a simple proportionality between the 
estimated and the 'actual' mutation rate is probably as transparent a relation 
as one could wish for. 

We are now ready to state our main 



Theorem 1 Suppose that the PGF for the number of mutants r in a large 
culture of bacteria is given as in £JJ). Suppose furthermore that the life-time 
distributions of neither mutant nor non-mutant bacteria are lattice-like, and 
that the mutation under consideration is neutral in the sense that the Malthusian 
parameter for a population of exclusively non-mutant bacteria is the same as that 
for a population of exclusively mutant ones. Then, with n\ as given in Equation 
fffi. there exists S E (—00, 00) such that, as the expected number m of mutations 
in the population tends to infinity, 

t r 1 . 5 
£m := ~i logy^m - 

converges in distribution to a random variable £ with characteristic function 
E(e^) = {-i9)- w iff condition (*) is fulfilled. 

We shall come to the proof in a moment, but let us remark now that the gist of 
the theorem lies in the scaling and centering of the distribution of mutants by 
n\m and lognjm, respectively, and in that the limiting distribution has been 
identified. This allows one to employ a maximum-likelihood method for esti- 
mating n\m and S/n' as described in [T2|. As of yet, the exact value of 5/n± 
does not seem to be of much importance, but the fact that, ultimately, the ex- 
periment will yield an estimate for the product of n* times the mutation rate 
rather than the mutation rate itself, does. We defer further discussion of this 
topic to the final section of the paper, and proceed to the 

Proof of Theorem^ By Lemma^J whose proof will be given in Section[S]below, 

m(g{s) - 1) =: m(l - s) {n\ log(l - s) + 5{s)) (10) 

for all s E [0, 1], where S(s) is bounded and has a finite limit as s — > 1. Therefore 
we obtain for the characteristic function ip m (6) of the random variable £ m , 

logVm(0) = m(l - e^K™)- 1 ) x 

(n\ log (1 - e^K™) -1 ) + 5(e tein ' mr1 )] - iO log n\m -i6-^ 

= -id log ( - i6(n\m)- 1 + 0(m- 2 )) + id -^ '- 

n' 

c 

—i9\ogn\m — i6 — + 0(m~ 2 ) , 
n i 

so that letting m —* 00, we find "0m(^) — * ip(Q) '■= {—i8)~ l6 if we choose 
6 := linis^! 5(s). This already concludes the proof of the theorem. □ 

Of course, the largest part work with the proof will be to show that the 
claim concerning the continuity of S(s) at 1~ made in (|10|l holds. The problem 
is to find a workable method of how to compare the function g(s) — 1 with what 



7 



is written on the right-hand side of (|TU)l . As we shall sec, there is one more 
or less natural way to do this if both mutant and non-mutant cells grow by 
way of an ordinary Galton- Watson process, but this is an example which does 
not meet the requirements of the theorem (the life-time distributions are very 
lattice-like in this instance), and we shall produce an example when a statement 
comparable to Theorem Q] does not survive in the Galton- Watson world. One 
must therefore exercise some caution so not to prove what isn't true. 

3 The Galton- Watson Scenario 

Suppose the generating function g(s) entering into Equation Q) were given by 

S(s):=: 9 f(s):=>^± f -^, (11) 

where fi(-) is for the z-th iterate of the generating function (@J, and fo(s) := s. 
We now prove the following 

Proposition 1 Let /(s) = s M for some integer fi > 2, and let g(s) be given as 
in Equation Ullp . Define 

7(s) ■-■Jf{s) := j— ■ 

Then lun s ^ 1 - "f(s) + Klog(l — s) does not exist in (— oo, oo) for any choice of k. 
Proof. One easily checks that 

W(s) = go f(s) + (n- l)s 

and 

VI (s) = : =n-l + h{s) 7 o f(s) . 

1 — s 

Let us now introduce a function $(•) and a change of variables s —: e~ e such 
that 

Thus, if f(s) = s^, we find 

0(t) = — + —*(*" log**), 

[i /i 

where hit) := h(e~ e ). Thus we see that over a range of log /x, i?(i) increases by 
about 1 — fi, so it makes sense to ask whether -d(t) might not be approximately 
linear in t; or, to be more precise, to ask whether 

r(t) := 0(t) - 4^ * = t(* " l0 S M) - — — 7(* - log /*) (12) 



has a limit as t — > oo. Because h(-) meets (*) if f(s) — s M , we have 

which implies that, indeed, limt^oo r(t) exists at least if t runs through integer 
multiples of log/x. This means that the lemma has been proved except for 
K = -t— — and if 

6(a) := 7(.s) + — log(-logs) 

had a finite limit as s — > 1~. We show that this is not the case: In fact, one 
easily checks that the coefficient of s k in j(a) equals fi~ l for /i i_1 < k < fi % — 1, 
and since /j, — 1 > log/x for any non-negative /z, it follows that 

. ,_ u 1W , M* Hlogfik 
- 1 5 U - lo S 



/zlog^ V ^log/"/ 

so that the sequence of partial sums of the first n coefficients of 6(s) cannot 
converge. However, these coefficients being of the order 1/fc for the fc-th coeffi- 
cient, they would have to if d(s) converged for s — > oo, by a Tauberian theorem 
of Hardy and Littlewood |3S] . But since 

log(- logs) - log(l - s) = log (^A -» 

as s — > 1, this contradiction concludes the proof of the lemma. □ 

In contrast to this, we have 
Proposition 2 Let 

m = 



s + /i(l - a) 
for some real number fi > 1, and let j(s) :=: 7/(s) be as in Proposition^] Then 

lim = -f(s) H logfl — s) 

exists in (— oo, oo). 
Proof. First observe that 

(/! - 1),U 



// — h(l — s) = s 



1 + (/J, - 1)8 



so condition (■*•) is fulfilled in this case, and that /(•) is of the fractional rational 
type. Hence we have for its i-th iterate |29| . 

1 - /<(*) 1 - 8 



fl l [L l + (1 — fl l )s 

which implies 

-. oo _, OO OO -, / -, \ fc 

Consider now the sum of the first n—1 coefficients. Except for the factor 1— /z , 
it is 

l / i \ r n-(i-/jT 



By Euler's formula, this equals 

f; i - (i - m~t = /°° (i - a - m~t) d« 

/•OO 

-nlog/i / (u- [u])/* _u (l - /i-")"- 1 du + 1 . 
./o 

We re- write the middle term as 

f°° ( nu,- u \ n ' 1 

log jj, / (u— [u])nu "1 dw 

Jo V " / 

J— log n/ log /j \ n / 

which readily follows after a change of variables u =: v + log n/ log /x. But this 
last expression tends to 

/oo 
(u — [u])^ - " exp(— ix - ") dv < oo 
-oo 

as n — > oo. Now one easily checks by expanding the integrand that 

/■CO IT 

/ (l_(l_ AI -)n) du= ^l. i 



where if„ is the n-th harmonic number. It follows that the discounted partial 
sums 



m i-i - 



H n 



fe =o,=o^V ^V lo ^ 
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converge, which by Abel's theorem already implies the lemma. □ 

This is something of a mess, because, as we have seen during the proof of 
Theorem^ the theorem actually asserts that the function 7(s) behaves essen- 
tially as — n* log(l — s) in the neighbourhood of 1, and by now we have both 
an example where it does, as well as an expample where it doesn't. The next 
step is therefore to suppose that g(s) is again as given in Equation (|TT|) . but 
with /(•) a PGF which isn't either s M or rrz—s • Then, with 7(5) as defined 
in Proposition Q] and with k = n~ , one finds that 

7 W + k io g (i - ,) = !tzl f; 1 i-^M (13) 

^ i=0 ' 

« v^ hofi(s) 1 - fi(s) ( , ,x 

"-g^? r -^log(/ l o/ i(s )) 

it z -^ ii l 1 — s v ' 

^ 1=0 ^ 

The way to see this is to use the definition J7J of h(-) and to check that the 
second and third sum, taken together, telescope. Now fix s 6 (q, 1), and denote 
by f-i(s) the inverse of f(s). Then /_„(s) := /_i o /i_ n (s) — » 1 as n — > 00 
(because /'(l) = /U > 1), and one deduces from lfT3|l that 

n 

--2 M i (l-/-i(*))fto/-i(*)log(fto/_ i (a)) 

n 

+ -J2n i (l-f- i (8))(^-hof_ i ( 8 ))]og(l-f- i (a)) 

^ i=i 

It can be shown [2U that under condition (*), 

/**(l -/_„(«)) -Q*(s) (15) 

asrn 00, with < Q*(s) < 00 for all s G (q, 1). Therefore, since 

A* 

the final right-hand term in the above equation is of the order 

C(i) 



E 
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and because (s — 1) log(l — s) = s + 0(s 2 ) for s in the neighbourhood of zero, 
we also have 

n n 

^2 [log [i - log (h o f-i(s))) + ^logho f_i( s )(n- ho f_i(s)) 



= i 



n _ n 

£ ( - l0 S C 1 - f-i(s))y a £( log (1 - f-i(s))) ~ £ 



C(i) 



Overall, we find that lining oo 7 o /_„(s) + ft log (l — /_„(s)j is of the order 
~ Y^i=i £{i)^~ lu) i 1 ~ a i which is bounded because of (•). The problem is that 
this only implies the existence of 

lim 7 o f_ n ( s ) + 4— L log (1 - f_ n (s)) 

n— >oo /Xlog/i 

for every fixed value of s 6 (q, 1), but doesn't say anything about whether this 
limit is independent of s or not. It turns out, however, that for Bellman-Harris 
processes with a non-lattice life-time distribution of particles, one can turn the 
above argument into an honest proof of Theorem Q] We will treat the example 
which is closest at hand in the following section. 

4 Markov Branching in Continuous Time 

We will now show that the argument given towards the end of the previous 
section goes through for continuous-time Markov models of clonal expansion. 
The construction goes as follows: Define 5' :=: 5/(3 for some 5 > 0, and let 

0(a) :=: <f>(s,S) := ^ J e^ u F u {s)du, (16) 

and 

( \ 1 x\ ^K 1 ) -4> oF S'{s) ,.„. 

V(s) ■■=■■ <PM ■■= (i) _fl 8 ) - < 17 ) 

Except for a norming factor, <j)(s) is a generating function, and we will find that 
<p(s) corresponds to the function h(s) defined in Equation (JJJ in a way which is 
sufficient for our purposes. Also, 

0(l ) = i(l-e^') = i(l-e-), 

4>'{1) = 1, and 



¥>(1) = lim - [ e -Pu l-Fu+6'(s) = 1 /■ e - f 3u e P { u+8>) du = e f38> = 

»-i 6' J 1 - s 6' Jo 
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provided that the Malthusian parameters for mutant and non-mutant cells are 
the same. This is where the assumption of no selection on mutant cells enters 
the calculations. Observe also that 

l-F u+s ,(s) = l-F 5 ,o F u (s) 5 
l-F u (s) 1-F u (s) 

so that e 5 — ip(s) > 0. Now, since the generating functions F u (s) have the 
semi-group property, we find 

/•oo /*<5 />oc 

(3 e- 0u F u {s)du = 13 e- pu F u {s)du + B e - f3u F u (s)du 
Jo Jo Js' 

= f35'cf>( S ) + -^ / e^ u F u oF 5 ,(s)du, 

eM Jo 

so that 

i=0 

where we have denoted by Fi*$(s) the i-th iterate of the PGF Fg(s). Note that 

^0(1) ,50(1) 

i=0 

as it must be. Also, 

l-g(s) _ ^ 1 ^-fFtfM 
7l ' Sj ~ 1-s -°2^ e iS 1 _ 8 

i=0 

and therefore 

7W + ^t^ ^ (w - *w) = *z 4 0(1) 'f:f" ,(s) d9) 

i=o e 

1 ^ y, o F itS ' (s) <f>(l) - o F^ (s) , . .v 

" ^ L ^ 731 log ^ ° F " 5 ' (s) ) 



gd Z_^ e ld l _ g 

4=0 

+ ^ 2^ ^ 1 - 1 log W 1 ) " ^ ° F « 5 ' (*)) • 

4=0 

Again, one checks this via inserting for log ip(-) into the second sum according 
to the definition of </?(•) and rearranging terms (which causes no trouble because 
of the absolute convergence of the sums), whereupon one ends up with a single 
sum which telescopes to 

0(l)-0oF o (s) , ■. 0(1) -<t>{s) i , , , 

log (0(1) - o F (s)) = - log (0(1) - 0(s)) . 



1-S ° V V ' " 1-8 
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The equality (|19fl established, it follows easily from there that, if we denote by 
F—i*S'(') the inverse of the function Fi*$>(-), 

7 o F_ n * 5 / (s) + i—= — log (0(1) - o F_„^/ (s)) 

1 ±* — n*5'\8) 

l-s 



e nS {l-F- n *s 
1 



( 7 ( S )+l O g(0(l)-0( S ))) 



£ e* 5 (0(1) - o F_*a,(*)) Ai(s) , (20) 



e nS (l-F_ n * 5 ,{s))j^ 
where 

A 4 ( S ) = 6 ^ oF J s '^ log (y, o F-i.,00) 

+ e '~^°J'" 3,(5) log (0(1) - o F_„„ («)) . (21) 

Thus, if we could show that the limit n — > oo in (|20|l exists and is independent 
of 5 and s, we would have made a big step towards a proof of the following 

Lemma 1 Under condition (*), and provided that the life-time distribution 
G'(t) of mutant cells is exponential (that is, G'(i) = 1 — e~ A * for some A > 0) 
we have that 

lim 7(s) + log(l — s) 

s— >1~ 

exists and is bounded away from ±oo. 

Proof. We first prove that the limit n — * oo in (|20|l exists. Because 



e 



(0(1) - o F_ w , (a)) = e 15 (1 - F_„ 5 , («)) ^ JY ^- 



converges as i — > oo (see [35], and observe that 0'(1) = 1), what we need to 
show is that the Aj(s) tend to zero sufficiently quickly. To do this, write M t for 
the expected size of a mutant clone of age t, and set H t (s) := ~._ ■ Then 
Af t = H t (l), and 

M t -ff t ( S ) 

= A/x / M t _ u e- A " du-xf ho F t - u {s)H t - u {s) e~ x ^ du 



Jo Jo 

= \ f (v~hoF t _ u {s))H t _ u {s)e~ Xu du 
Jo 

+ A/i / (Af t _ u - #«_«(«)) e" Au du , 
Jo 

which can be solved by means of Laplace transforms: 

M t - F 4 (s) = A / (n-ho F t _ u (s))H t _ u (s) e 0u du . 
Jo 



14 



Then, since M t = e 131 , it is straightforward to check that 

< e 5 - ip(s) < i J e-P u (M u+s , - H u+S , («)) du 

= J, j e ~ Pu I (n-hoF tt+s >- v (s))H u+sl - v (s)e^dvdu 
< \(fj,-hoF 2S ,(s)). 

But it follows easily from Jensen's inequality that F t (s) > s Mt — s e , which by 
(*) immediately gives 

e 5 - <p{l - s) < \e 2Su) s"( - logs - 2<5)~ Q £( - logs - 26) . 

This is good enough for us: because 

1 ( j? M \ j,, f, e s - (p o F-„ s > (s) \ 
log (ip o F-„ s > (s)) = 6 + log I 1 - s I 

differs, in the neighbourhood of 1, from 5 only by an amount of order e — 
(f o F_i»5'(s), the Ai(s) as given by (|21|l tend to zero at least as quickly as 
~ C(i)i l ~ a , which proves that the sum in Equation i|20|) converges. But 



f 



^^ = U S 'e-^(M u -H u ( s) )du. 

1 — s o Jo 



which by the above implies that 

<Ki)-<Ks) 



log(l - s) - log(l -s)->0 
1 — s 

as s — > 1. All in all, we know now that 

lim 7 o F^ n *s> (s) + log (1 - F_ n * s > (a)) (22) 

n — >oo 

converges for every <5' > and every s £ (9,1)- But i*t(-) is continuous in 
i; and since 'y(s) and log(l — s) are continuous for s G [0, 1), an application of 
Kingman's lemma [3S] shows that the limit in lll'l'l) is independent of 8. Since we 
can find for any s < u, s,u £ (q, 1) some suitable n and 5 such that s = F n *s> (u), 
this limit is also independent of s, which concludes the proof of Lemma ^ D 

4.1 A Second Proof 

For completeness, we now give a second proof of Lemma ^ which makes use of 
the forward picture for the continuous time Markov process, namely 29 

|*(-) =«(-)£*(-), 
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where 

u(s) = \{f(s) - s) . 

Then we have by definition of g(s), 

oo poo o 

g(a) = -e-^F u {s) + e -^—F u {s)du 
u=0 J ou 

f°° 8 X 

= S -X(s-f(s))j e-^—F u ( S )du = s--( S -f( S ))g'( S ), 

which, as one readily checks, implies 

7(*) = l+5(fl)(7(*)-V (*)(!-«)) (23) 

where 

p 1 — S fJ, — 1 

Now define t :— — log(l — s) and a function #(•) of t such that j(s) —: #( — 
log(l-s)). Then 

0(t) = i + *(*)(<>(*) -0'(t)), 

and z(i) := z(l — e — *). This can be solved by standard methods. We have 

1 



where 

Mt) ■■= m^V (f ^h 1 du 

is the solution of the homogeneous equation corresponding to ill'.'il) . and q G 
(q, 1). Now write 

J s ■d {u)z(u) J,. z(u) 

and observe that $o(t) — #o(w) is at worst of order 

/"' z(V) - 1 
0><?o(*)-i?o(tt)> / , x d« 

y« z ( w ) 



in magnitude, so 

/' 



>f!|HoM^ _i_r (u _ £(!0^ du . 
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which by definition l|24|) of z(-) and because of (*) is of order 

- [ (u- ^/7 f , du ~ - /" e- uu C(u)u a - x du > -oo . 

7, v ; ft(i - e-«) - 1 y e w 

This concludes the proof, since now it follows that 

7(a) = 0( - log(l -*))=- log(l - a) + <5( - log(l - «)) , 
and <5(-) approaches some finite value as its argument tends to infinity. D 

By now, we have completed the proof of Theorem ^ m case the life-time 
distribution of mutant cells is exponential. From a practical point of view, this 
is what we wanted the least, since in a real-life situation, newborn cells will have 
to pass through a whole cell cycle before they can divide anew. Therefore, in 
the following section, we shall imitate the first proof of Lemma ^ to obtain a 
comparable statement for general models of cell proliferation. 

5 The General Case 

One difficulty with non-exponential life-time distributions is that, if one only 
looks at the number of particles Z t at each given instant, the process {Z t }t>o 
is not Markovian anymore. This implies that the generating functions F t (s) = 
E(s Zt ) will not form a semi-group, so that in particular, F_ u+t o F_ t (s) ^ 
F_„(s). Nevertheless, limt-xx, F^ u+t o F_ t (s) exists and can be calculated as 
follows: Fix s 6 (q, 1), where g, again, is the unique fixed point of f(s) in 
[0, 1). It can be shown [3U| that the product of — log F_ t (s) times the number 
of bacteria Z t that by time t have originated from a bacterium born at time 
converges almost surely to a non-degenerate random variable Z. Then, with 

R(x) :=: RJx) := lim F t o F_ t (s) x (25) 

t— i-OO 

the Laplace transform of the random variable Z, we have [50] 

/•OO 

lim F_ u+t o F_ t {s) x = Rixe-^) = / o i?^^) dG(u) . (26) 

t->ao J 

Consider now the number of bacteria Z t [x;y) aged less than y at time t that 
have originated from a single bacterium aged x at time 0. Then 



Z t [x]:= / Z t [x;dy] (27) 

Jo 

is for the total number of progeny an individual aged x has produced by time t, 
and it is easy to see (say, by considering the corresponding Laplace transforms) 
that — log .Fit (s) times this number converges in distribution to random variable 
Z[x\. Now for r(x) a measurable function on [0, oo), define 

Zt[x] 



/■00 J- 1~ J 

/ r(y)Z t [x;dy] := V r(y 

•>0 i-i 
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where the yi are the ages of individual bacteria in the population at time t, and 
set 



§t[x](r) :=Ef exp(f™logr(y)Z t [x;dy]) 1 . (28) 

This is the definition of the generating functional (GF). The 3 > t[-](-)' s do have 
the semi-group property in the sense that <&t+ r [x] (r) = <& t [x](<Jv[-](r)). Fur- 
thermore, Schuh |30j has shown that there exist measurable functions s t (x) on 
[0, oo) such that 

$ r [x}{s t )=s t _ r (x) (29) 

for all < r < t < oo. The S{(x)'s can be specified explicitly [30]: With 

the residual life-time distribution function of a mutant bacterium aged y (we 
write G(i) instead of G'(t) for brevity), 

/•OO 

s t (x):= foR( e -^ t+ ^)dG x (u). (31) 

Jo 

Since R(0) = 1, the s t (x) tend to 1 as t — > oo. We need an estimate for the rate 
of convergence. From now on, we will always assume that G(-) is non-lattice. 
We have the following 

Lemma 2 Fors t (x) as given by Equation ijffij) . -R(-) as given in Equation $25\) . 
and 

/■OO 

V(x):= / e- 0u dG x {u) 
Jo 

the reproductive value, we have 



for some constant C > and any u < t. 



<^-hoR( e -^ t - u '>) 



The proof of this lemma is deferred to the appendix. Observe that convergence 
is uniform in x. Next, we prove 



Lemma 3 For 



POO 

g[x](r):=0 e- f3u $ u [x}(r)du, 
Jo 



with <fr u [x}(r) as given in Equation \2tft) . and s t (x) :— exp (n* log Sj (x)j , we 
have that 

lim _ gwfl^frjw _ , r (e _,„ M4x] _ „.„ nx)) du . 



logi?(e"' 3 *) 
In particular, the limit exists in (—00,00). 
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Proof. Note first that the lemma is trivial if Z u [x] explodes in finite time. We 
may therefore assume that Z u [x] < oo for all t £ [0, oo). Also 



limZr 1 / Z t [0;dy]-+A(y) 



almost surely, where A(-) is the limiting age distribution of particles 30 . Granted 
this fact, the same is obviously true for Z t \x\~ x JJ? Z t [x; dy\. Now for the differ- 
ence u — t fixed, we have 
i 
logs t (y)Z u [x;dy] 

log s t (y) logi?(e-' 3t ) Z u [x;dy] 



= nllogRie'^Z^x] 



lo gj R(e-^) logi?(e-' 3 «) Z u [x] 



-► - n?ixe^ u -^Z\x\ / V(y) dA(y) = - n\e^ u -^ Z[x\ 
Jo 

at least in distribution as t — > oo, and one finds 

t^ooj t logi?(e-P*j 

V ; t-Wt logi?(e-^*) 

Now e _a: is a completely monotonous function in x; thus e~ y (y — x) < e~ x — 
e~ y < e~ x (y—x) for arbitrary a; and y. Because of this, and because of Equation 
l|32|). we can say that 



,. G[x]{R(e-^))-g[x](s t ) 

lllTl : ; 3TT 

*^oo logi^e - '") 



log s t (y) 



= lim /3 / e"^ E £ t (u)Z u [x] + n\£ t (u) / L _i» Z u [x; dy] du , 

where Z u [x] is as in Equation (|27() . Now £t(u) is of order exp(— e~^*Z u [a;]). 
Thus, £t{u) tends to 1 as t — > oo for every fixed value of u, and the lemma will 
be proved once we have shown that 

13 J* e-^Ef Z u [x] + n{ J°° J°^L ^.[*i M J <*« 

converges as £ — > oo. By linearity of taking expectations, the above can be 
re- written as 

lim f3 / (e~ fiu MJx] - nlfjV(x)) du 

'^°° Jo ^ ' 



lu=0 Jy=0 

19 



since V{y) is a right eigenfunction to the kernel M t [x;y] := E(Z 4 [a;;y]) with 
eigenvalue e' 3 ' (3UJ. But 

e -fr* t -<*£(fi f e - f3u M u [x]du 



e 

u=0 Jy=0 



because of Lemma [21 and this tends to zero for large values of t. It therefore 
remains to show that e~ l3u M u [x] — n\^V{x) is integrable. We have [30] 



( 



-m 



M t [y}-nl t iV(y) = e- f3t (l-G y (t)) 



e~ 0u dG v (u) + /i / (e- 0{t - u) M t _ u -n\)e- 0u dG v {u) 



n x n 



where Mt is short for Mt[0], and one readily checks that 



OC 



e 



-fit. 



l-Gy(t))dt 



1 - V(y) 



/o P 

and 



OO /«OC 



e-P" dG v («) dt = / w^" dG y (u) < 

t=QJu=t JO 



Furthermore, 
ft 



t=o J«=o 



(• e -/3(f-«) Mt _ u _ n .) e -/9u dGy ( w ) ^ = v(y) / (e-^Mt - nj) dt 



so that the integrability of e ^ t Mt[y]—n'fj,V(y) follows from that of e /3t M t — n°. 
But the latter is just Proposition 8 in [37]. □ 

Lemma|21fails in the Galton- Watson scenario, and the reason for this is just 
the lack of integrability of e~^ l M t — n* in this case (n* should be taken equal to 
the value of k in Equation (JT3JI). This is not in contradiction with Proposition 
[21 it only shows that with the methods we use for the proof of, ultimately, 
Theorem ^ we cannot prove or disprove either of Propositions ^ an d 

We now come to the main result of this section: 

Lemma 4 Suppose that the life-time distribution G* (t) of mutant cells is not 
lattice-like, and that g(s) is as given in Equation 0) with lo(s) — s. Then, with 
n* as defined in Equation |^). we have that 

1 — q(s) 
lim - yy ' + n\ log(l - s) 
s->l- 1 — s 

exists and is bounded away from ±oo. 
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Proof. First observe that g(s) = Q[0](s), where s(x) is that function on [0, oo) 
which assumes the constant value s for every x, and £/[0](s) is as given in Lemma 
13 Therefore it will be sufficient to show that 

m^+"i*<-**c-*» 

converges as t — > oo. By Lemma [21 this boils down to proving convergence for 



logi?(e-/ 3 ' 



log(-logi?(e-^)), (33) 



and to demonstrate this, we define 

^](r):=:/j[i](r):=i/ e-" u * tt [a;](r)du (34) 

and 

WWW :=: ««NW-^pi^- • ( 35 ) 

We need to know about the behaviour of these two functionals in some neigh- 
bourhood of the unit function 1. First, JF[x](l) = (1 — e~ s )/S. Now choose any 
t > and observe that, because of Equation l|29|l . 

logfl(e-/") = W„ 6 logfl(e-^) d " " "^ (36) 

as £ — > oo. Moreover, by Lemma El the difference between the two sides is 
bounded in absolute value by some constant times \i — ho R(e~^^~ s )). Next, 



f!L f S \-0(u+s>) l-st-u-8>(x) logRje-*) 



fin 
5' Jo ' logi?(e-^) ^](l)-^](s t ) 



e s 



iff mutation is neutral, and convergence is at least as quick as ~ [/, — h o 
i?(e- /3( *- 2,5 ' ) ), again by Lemma[2 

Returning now to the proof of Lemma QJ it is readily checked that, by the 
definition of the functional C/[ir](r) in Lemma|31 

Q[x]{v) = 8^ — s . (38) 

i=0 

Thus, we obtain in complete analogy with Equation (|19|) , 
1 -£?[*] (r) , ^](l)-^](r) 



logi^e"/ 3 *) logi?(e-^) 
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log(^](l)-^](r)) (39) 



" 1 ^[x](1)-J=[x](^ s ,[.](t)) 
2-s e iS logRie-Pt) 

1 » / «H(^H(r)) J=[x}{l)-Hx](^M-](r)) 
e& ho\ e ' 5 logi?(e"' 3 *) 

log («[!](*(,,/ [-](r)) 

1 ~ / ^-HH^Hr) ) ^](l)-^](4W[-](r)) 

/ ■■ I P iS 



e s 



i=0 



e lS logRie-P*) 

log(^c](l)-^B](^[-](r)))V 

By a similar line of argument as for the proof of Lemma ^ we find that the 
expression on the right-hand side of the above equation, when taken at r(-) = 
s t (-), converges as t — » oo through integer multiples of 5' . As an illustration, 
consider the final right-hand term in the above array, for which we claim that 

1 y,/ V-tt[a:]($i.y[-](B )) J=[x]{l) - J=[x](^s>[-Kao)) 
e sA^\ e i8 e nS \ogR(e~ nS ) 



log(^](l)-^]($«y[-](s )) 



1 - 

7s Z^ 



i<5 I 



„$ 



W[a;](si* g /) ^[»](1) - ^t*](si*5/) 



e S Z_^ l e n5 logi?(e-™' 5 ) 



j=l 



]0g(f[x](l)-f[x](Bi 



converges as n — > oo. But this is clear from Equations 136J I. 137(1 . assumption 
(*), and the fact that e l3t logR(e~ /3t ) approaches a finite limit as t tends to 
infinity 30 . 

Thus we have shown that 



n — >oo 



\ogR(e-^) + logR{e-"*) M"^^ > 



Fs[x]{l) - Fs[x](s n s>) / Jg[x](l) - .Fg [x] (s n g/ ) 
logR(e- nS ) ° g \ -logR(e- nS ) 

exists for fixed S. We now have to deal with the fact that .^['K') depends on 6 as 
an external parameter. Because of (1361) . the final term in the above expression 
has a finite limit which does not depend on <5, and from the little discussion 
following 136|l it is clear that for < r < S the difference 

■fr[g](l) - -FrN(St) fs[x](l) - Ts[x](s t ) 

log R{e-I 3t ) \ogR{e-' 3t ) 
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can be of order at worst /j, — h o R[e /3( - t 5 - ) ). It follows that 
also exists for fixed r and arbitrary S > Now 



and 



G[x](at) = /3 / e-^s t _ u (aO dw + e^'/? / e~' 3 "$ u [x](s ) cfc 
Jo Jo 



R p -Pt ft 

T T [x]{ St ) = l — / e Pu s u du 

T Jt-r/p 



are continuous in t, and the same is obviously true for R(e "*). Thus it follows 
again by Kingman's lemma that 

t^oo log R(e p 1 ) log it(e P 1 ) v ' 

exists in (—00, 00), which by Equation l|36|) implies that the same is true for 



l-G[x](B t ) 

t-T<i> \ogR( e -' 3t ) 



!™ L5 - ^W lo S ( " lo g^( e ^)) 



But this proves the claim made for J23J, since /j,V(0) = 1, and thus the proof of 
Lemma0]and, finally, that of Theorem^is complete. □ 

6 Discussion 

The good news to take home from this paper is that, under the condition that 
m is large enough for most cultures in a fluctuation experiment to produce a 
fairly large number of mutant colonics (by comparison with Lea and Coulson's 
[H?| results, one would venture that to > 4 is already sufficient), it is possible 
to forget about the subtleties of clonal proliferation and still come up with an 
estimate for the mutation rate that makes better use of the data than just 
looking at the fraction of cultures without any mutants. Indeed, starting from 
Lea and Coulson's observation that for suitably chosen constants a, b, and c, 
the variate a/(£ + b) — c (with £ as in Theorem ^) is approximately normal 
distributed, one would go about evaluating the experiment as exemplified in 
|T§|. or rather (because of the additional summand 8/n') as developed further 
in m>| . The advantage of this method is that additional factors like phenotypic 
lag of mutants, or residual growth of non-mutant bacteria on selection medium, 
readily integrate into the calculations. The bad news is that the estimate for m 
is consistently wrong by a factor n' . This confirms an observation of Oprea and 
Kepler jS], who have found that their continuum approximation of the Luria- 
Delbriick distribution can be adapted to realistic cell-cycle time distributions by 
a two-parameter generalisation. Although Oprea and Kepler are only concerned 
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with the case n = 2, one can see that their factor 6/2 corresponds to our n', and 
it seems that the natural logarithm of their factor c corresponds to our S/n'. 

It is natural to ask about the values of n* and 6 for realistic models of cell 
proliferation. Since S/n* can in principle be estimated from the experimental 
data m^! the question of practical relevance really is about the magnitude of 
n*. For Kendall's ^H] multi-stage model of cell proliferation, for instance, with 
the bacteria lingering in each of k stages for an exponentially distributed length 
of time, the life-time distributions are Gamma-like with densities 



k \k 



k k X 



kXtj.k-1 



r(fc) 

which implies (3 — k\(n k — 1), and 



e~ KA H 



1 



k(fi — fi 1 k ) 

so that n\ approaches 1/fc for very large values of /i. For Rahn's [HH! life-time 
distribution density 

ake- at {l-e- at ) k - 1 , 

which obtains under the assumption that cells linger in each of in total k stages 
for a time that is Gamma distributed with shape parameter 2, the respective 
parameters cannot be evaluated explicitly anymore (at least not unless k < 4), 
but it seems that (3 « a(k + l)(^ fe — l)/2, and 

1 / k 1 v _1 

n, = — 



E 



/3 \x \~[ (3/a + i I 

It is not difficult to prove via Holder's inequality that for every life-time distri- 
bution of bacteria, 

n i > —, > 

fxlogfi 

but (by taking for G'(-) the Gamma distributions) it is also not difficult to see 
that n* cannot be bounded from above. Still, it would seem that for realistic 
models of clonal proliferation, n* varies in the range from 0.7 to 1. This would 
result in an underestimation of the mutation rate by, at worst, 70 percent, which 
isn't very much, but for experiments with a sufficiently large number of cultures 
(say, about 100) could differ significantly from an unbiased estimate. 

The question about the magnitude of S seems interesting in itself. The results 
of Oprea and Kepler suggest that at least for the cell-cycle time distributions 
investigated in their paper [SJ, 8 is positive and of about the same magnitude 
as n*. The proofs of Lemmas ^ an d E partially confirm this. In fact, since in 
practice the offspring distribution of (mutant and non-mutant) bacteria has mo- 
ments of any order, the u> in condition (*) can be taken as positive, so that the 
sums in Equations l|19|l and (|39(l converge geometrically, and are about 1 — /i _1 
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in magnitude. Finally, it seems possible to demonstrate some 'stability' of the 
Luria-Delbruck distribution with respect to different models of clonal prolifera- 
tion also if the assumption of neutrality of the mutation is dropped, although, 
unsurprisingly, one ends up with stable limit laws of index other than 1. It is 
planned to present these results in a future publication. 
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7 Appendix 

Our goal in this appendix is to supply a proof of Lemma |2 Because 
1 1 _ 11 + |(1 -a) + ... 



logs l-« 21 + i(l-s) + ... ' 

we may as well prove the corresponding statement for 1 — i?(e~^*) in the de- 
nominator instead of — logi?(e - ^*). We have 



fj,V(x) — e 



-0u 1 s t-u{X) 

1 - R(e-^) 



1 - R{e-^) 

/ { l x-hoR{ e -^ t - u+v ^))e' 0v dG x {v) 
Jo 

- r h o R(e-^-^) X{t ~ U + V J ~ X{t) e-»> dG x (v) , (40) 



where 



X{t):=e fit {l-R{e' 0t )), (41) 

and R(-) is the Laplace transform from Equation I|25J1. By monotonicity of h(-), 
the first term on the right-hand side of the above Equation (|40() is no larger 
than 

V(x)(n-hoR(e-W-^)) . 

The problem is hence to obtain an estimate for the second term which, as X(t) 
is increasing in t, consists itself of a positive and a negative part. It is however 
easy to see that neither of these can be larger than 

f" h o fl(e-^-"+»>) X{t - U t V) " X{t " U) e~ 0v dG x (v) , (42) 

Jo x (t ~ u ) 
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so we may as well attempt to bound this. Furthermore, since it is only the ratio 
X(t — u + v)X(t — u)^ 1 which enters Equation l|42[l . we may, for ease of writing, 
assume that X(t — u) — 1 for < t — u arbitrary, but fixed. We now deduce 
from Equation (|26|l that 

/>oo 

X{t) = hoR( e -^ t+u ^)X(t + u)e' fitt dG{u) ) (43) 

Jo 

and then, if we write 

G (t) :=n J e-P u dG(u), (44) 

Jo 



that 



X(t)= I k - fi(e U) — ^X{u)dG ti {u-t)+ I X{u)dGp{u-t) 



This is a relation which can be iterated: we use l|43[) to express the second 
integral in the above equation as 

x{t) _rhoR{e-^)-, 

it V 

I '-X{v) dG p {v - u) dG p (u - t) 

u—t \J v—u ft / 

h o R{e-P v ) 



X(v)d[ / G (v-u)dG (u-t) 

v=t A* \Ju=t 

^^ >-X(t + v) dGf(v) , 

o A* 

where G*g{t) is for the two-fold convolution of Gp(t) with itself, and proceed 
by induction: 

X(t)= / X(t + u)dG% n (u) 
Jo 

r°° a - h o R( e - 0< - t+u h -A 

- / ^ [ >-X(t + u) Y, dG% (u) (45) 

for every natural 1 < n < oo. Now the 'renewal function' 

oo 

Up(t):=J2 G t^ ( 46 ) 

(the standard definition would be to start the summation with i = 0) is finite 
on bounded intervals [2HJ, whence it follows that 

/•OO 

lim / X(t + u) dG* g n (u) = lim X(t + u) =: X x , 

n— >oo Jq " u-^-oo 
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where the limit exists in [0, oo] because the other terms in Equation l|45|) are 
obviously monotonous in n. Assumption (•) implies that, in fact, X^ < oo |30|. 
We finally obtain 

X oc -X(t) = / — '-X(t + u)dU (u). (47) 

Jo M 

By renewal theory J2] > 
for t > 0, where 

/■oc /-oo 

z^ 1 := / tdG p (t) =/j, te- f3t dG(t) , 
Jo Jo 



and U(-) is a measure which (is concentrated on the nonnegative reals and) 
integrates to 

i= & f°° t 2 e -,t dG(i) 

Thus 



z Jo 



#„, - X(t) = v 8 r M hoR ( e -^ Af (u) du 

+ / £ !s J -X(t + u)dU p (u), 

Jo V 

and because h o R(e~ l3t ) < /U for t 6 [0, oo), it then follows that 

ft o i?( e -^'-"+ 1 ') ) (# (t - u + v) - 1) e- f3v dG x (v) 

ii 

/■CO rt — U-\-V 

< vp I I (ix-ho R(e- 0W )) X{w) dv e~ 0v dG x (v) 

J v—0 J w—t — u 

+ / e-f* v dG x (v) ( f i-hoR(e~ fj ( t ~ u+w '>))X(t-u + w)dU l3 (w) 

Jv=0 Jw=0 



< 



X(0) 



(v-hoR(e-W-^))(apV(x) + Vp J ue^" dG x (u)\ , 



where we have divided by X(0) to get rid of the assumption X(t — u) = 1. This 
obviously concludes the proof of the lemma, since V(x) < 1, and the integral is 
bounded by (e/3) -1 . □ 
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